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1. Summary 


In April 1995, the University of Illinois received a NASA Grant NAG3-1474 entitled 
“Electromagnetic Scattering from Realistic Targets” for a one-year period and the grant was 
extended for another year with no cost. The co-principal investigators are Professors 
Shung-Wu Lee and Jianming Jin, and the NASA technical monitor is Dr. Afroz Zaman. 
This project has been completed. Results are reported in the following technical reports: 

1. J. M. Jin, J. A. Berrie, R. Kipp, and S. W. Lee, "Calculation of radiation patterns 
of microstrip antennas on cylindrical bodies of arbitrary cross section," May 1996. 

2. A. D. Greenwood, S. S. Ni, J. M. Jin, and S. W. Lee, "Hybrid FEM/SBR method 
to compute the radiation pattern from a microstrip patch antenna in a complex 
geometry," October 1996. 

3. G. Fan and J. M. Jin, "Scattering from a cylindrically conformal slotted-waveguide 
array antenna," December 1996. 

4. G. Fan and J. M. Jin, "Scattering from a large planar slotted waveguide array 
antenna," August 1997. 

5. E. Branch, "A numerical comparison of measured SAR images with XPATCH 
computed SAR images," UIUC EM Lab. Scientific Report No. 95-2, September 1995. 

6. S. Ni, J. M. Jin, and S. W. Lee, "Hybridization of the finite-element method and the 
shooting-and-bouncing-ray method for scattering and radiation from large and complex 
targets," UIUC EM Lab. Scientific Report No. 95-3, October 1995. 

The first four reports are included here as appendices. The last two reports have been 
submitted to NASA before and, hence, are not included here. 

The general goal of the project is to develop computational tools for calculating radar 
signature of realistic targets. During the period of this grant, we developed a hybrid 
technique that combines the shooting-and-bouncing-ray (SBR) method and the finite- 
element method (FEM) for the radiation characterization of microstrip patch antennas in a 
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complex geometry. The technical details are given in Appendices 1 and 2. In addition, we 
developed a hybridization procedure to combine the moment method (MoM) solution and 
the SBR method to treat the scattering of waveguide slot arrays on an aircraft. The technical 
details are described in Appendices 3 and 4. 


2. Publications Resulting from This Research 

Journal articles (published or accepted): 

1. J. M. Jin, J. A. Berrie, R. Kipp, and S. W. Lee, "Calculation of radiation patterns of 
microstrip antennas on cylindrical bodies of arbitrary cross section," IEEE Trans. Antennas 
Propagat., vol. AP-45, no. 1, pp. 126-132, Jan. 1997. 

2. A. D. Greenwood, S. S. Ni, J. M. Jin, and S. W. Lee, "Hybrid FEM/SBR method to 
compute the radiation pattern from a microstrip patch antenna in a complex geometry," 
Microwave Opt. Tech. Lett., vol. 13, no. 2, pp. 84-87, Oct. 1996. 

3. G. Fan and J. M. Jin, "Scattering from a cylindrically conformal slotted-waveguide 
array antenna," IEEE Trans. Antennas Propagat., vol. AP-45, no. 7, pp. 1 150-1 159, July 
1997. 

4. A. D. Greenwood and J. M. Jin, "Hybrid MoM/SBR method to compute scattering 
from a slot array in a complex geometry," Applied Computational Electromagnetics Society 
Journal , accepted for publication, 1997. 

5. G. Fan and J. M. Jin, "Scattering from a large planar slotted waveguide array antenna," 
Electromagnetics, accepted for publication, 1997. 

Conference papers: 

1. G. Fan and J. M. Jin, "Scattering from a cylindrically conformal slotted-waveguide 
array antenna," The 1996 IEEE Antennas and Propagation Society International 
Symposium, Baltimore, MD, 1996, pp. 1394-1397. 
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2. A. D. Greenwood, S. Ni, and J. M. Jin, "Computation of the radiation pattern of a 
microstrip patch antenna in a complex geometry," The 1996 IEEE Antennas and 
Propagation Society International Symposium, Baltimore, MD, 1996, pp. 256-259. 

3. A. D. Greenwood and J. M. Jin, "Hybrid MoM/SBR method to compute scattering 
from a slot array antenna in a complex geometry," The 13th Annual Review of Progress in 
Applied Computational Electromag-netics, Monterey, CA, 1997. 
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APPENDIX 1 


Calculation of Radiation Patterns of Microstrip Antennas on 
Cylindrical Bodies of Arbitrary Cross Section 


Jian-Ming Jin, Jeffery A. Berrie, Robert Kipp and Shung-Wu Lee 


ABSTRACT 

A simplified approach based on reciprocity is presented to calculate the radiation patterns of 
microstrip patch antennas or arrays on a cylindrical body having an arbitrary cross section. 
In this approach, the microstrip patch antennas are characterized using the finite element 
method, and their radiation patterns are then calculated using a two-dimensional method of 
moments in conjunction with the reciprocity theorem. The validity of the calculation is 
demonstrated by comparison with measured data for a microstrip patch antenna on a 
circular cylinder attached to a plate. Other numerical examples are also presented to show 
the capability of the method, as well as various effects of the host cylinder on radiation 
patterns. 


This work was supported by the Office of Naval Research under the grant N00014-95-1-0848, NASA under 
the grant NAG3-1474, and the Air Force Rome Laboratory under contract F-30602-93-C-0036. 

J. M. Jin and S. W. Lee are with the Electromagnetics Laboratory, Department of Electrical and Computer 
Engineering, University of Illinois at Urbana-Champaign, Urbana, IL 61801-2991. 

J. A. Berrie is with Mission Research Corporation, 3975 Research Boulevard, Dayton, OH 45430. 

R. Kipp is with DEMACO, 100 Trade Center Drive, Champaign, IL 61820. 
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I. Introduction 


Microstrip patch antennas [l]-[3] are attractive for aircraft and spacecraft applications 
because they are conformal to the surface of the host object. However, their analysis is 
often performed by assuming that the antennas are placed on an infinite ground plane with 
an infinite dielectric substrate [4]-[6] or on a finite substrate housed in a cavity recessed in 
the ground plane [7]. Although such an analysis can be extended to the antennas on a 
coated circular cylinder [8]-[ 1 0] or in a cavity recessed in a conducting circular cylinder 
[11], it cannot be applied to those on a cylinder with an arbitrary cross section because the 
Green's function for such a structure is not available. While the shape of a host cylinder 
has little effect on the input impedance and current distribution of a microstrip antenna since 
the latter is a highly resonant structure, it can alter the radiation patterns substantially. The 
objective of this paper is to develop a simple method to calculate the radiation patterns of 
microstrip patch antennas on substrates residing in cavities recessed in cylindrical bodies of 
arbitrary cross section possibly with dielectric coatings. 

In the proposed method, the equivalence principle is first employed to divide the problem 
into two equivalent problems. The first problem is to characterize microstrip patch antennas 
housed in a cavity whose aperture is closed with a perfect conductor, which is 
accomplished using the finite element method [7]. The second problem is to compute the 
field radiated by the equivalent magnetic current over the cavity's aperture in the presence 
of the host cylinder, and to accomplish this task, an approach based on the reciprocity 
theorem is developed in conjunction with a two-dimensional method of moments for 
scattering computation. The validity of the proposed method is demonstrated by 
comparison with measured data for a microstrip patch antenna on a circular cylinder 
attached to a plate and for a waveguide-fed trihedral. Other numerical examples are also 
presented to show the capability of the method and the effects of patch rotation, material 
coating, and finiteness of a ground plane on the radiation patterns. 


II. Formulation 

Consider the problem illustrated in Fig. 1 where several microstrip patch antennas are 
placed on a conducting cylinder of arbitrary cross section. These microstrip patch antennas 
are housed in one or several shallow cavities and the conducting cylinder may be coated 
with dielectrics. Since a microstrip patch antenna is a highly resonant structure, its current 
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distribution and input impedance at resonance are mainly determined by internal structures 
such as the shape of the patch, the thickness and dielectric constant of the substrate and 
superstrate. Provided that a patch antenna is placed on a locally flat surface or surface with 
a small curvature, the shape of the host cylinder (external structure) has little effect on its 
current distribution and input impedance, as demonstrated in [11]. Therefore, the patch 
antenna can be characterized approximately by assuming that it is placed on an infinitely 
large ground plane. Applying the equivalence principle, we can close the aperture of the 
cavity housing the patch antennas with a sheet of perfect conductor and introduce a 
magnetic current density over the aperture defined as 


M = Exn 


( 1 ) 


where h denotes the outward unit vector normal to the aperture. As shown in [7], this 
magnetic current can be determined using the finite element method with vector elements by 
seeking the stationary point of the functional 


F(£) = -|jJ 

* \ 

+1U 


1 


(V x £) • (V x £) - k.le r E ■ E 


dV 


V L r~r 


1 .-.i 


jk 0 Z 0 J' nt ■ E - — M ,n ‘ • (V x E) 
f^r 


-* 0 2 jjM(F)- \\G 0 (r,T') M{r’)dS' 


dV 


dS 


( 2 ) 


where V denotes the volume of the cavity and S denotes the aperture of the cavity. Also, 
G 0 (r,r') denotes the free-space dyadic Green's function, and J ,nt and M ,m represent the 

internal sources or antenna excitation. The procedure of the finite element solution and 
modeling of a current probe and impedance load are described in [7] and [12], and the 
modeling of waveguide feeds is described in [13]. 

Once the equivalent magnetic current M over the aperture of the cavity is calculated, the 
radiation patterns can be obtained by computing the field radiated by M in the presence of 
the conducting cylinder with the cavity filled with perfect conductor. Since the cross section 
of the cylinder is arbitrary and the cylinder may be coated with dielectrics, no pertinent 
Green's function can be derived and the radiated field cannot be expressed by an integral 
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formulation. To alleviate this difficulty, we propose to employ the reciprocity theorem to 
calculate the radiated field. According to the reciprocity theorem [14], we have 

JJJ E a ■ J h dV = -JJJ H h ■ M“dV (3) 

where E“ is the field produced by M“ and, likewise, H h is the field produced by J h . 
Letting M a be the equivalent magnetic current M over the aperture of the cavity and J b be 
the electric dipole of moment II located at (/-,0,0) and oriented in the u -direction, the 
radiated field is then given by 

£(r,0,0) w = --JJtf h MdS (4) 

Often, the radiation patterns in the plane perpendicular to the cylinder (that is, the plane of 
6 = n/2) are of particular interest. In this case, we obtain 

E z (r, 0) = jkoZ / - f f H b M r dS (5) 

4 nr 

lie £ r r 

E t (r,<t>) = - j f H h z M z dS (6) 

47rr 

where H h r (or H z ) is the magnetic field induced on the surface of the cylinder (again, with 

the cavity filled with a perfect conductor) by an incident plane wave whose electric (or 
magnetic) field is in the z-direction and has a unit amplitude, and M T (or M z ) is the 

component of M in the direction perpendicular (or parallel) to the cylinder. 

It is then clear that to calculate the radiated field, we have to find H b and/or H z first. This 
can be accomplished simply by solving the problem of plane wave scattering by the 
cylinder using the method of moments. There are several different formulations available 
for this problem. The most efficient one is to employ the volume-surface integral equation 
[15] given by 
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for -polarization, and 
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for // z -polarization, where Q denotes the cross section of the dielectric coating, 
denotes the interface between two different dielectric materials, and r r denotes the 
conducting surface. Also, n d denotes the unit vector normal to r d and pointing from the 
" - " side to the " + " side, and G 0 represents the well-known two-dimensional free-space 

Green's function. For cylinders without coatings, (7) and (8) reduce to the well-known 
equations 
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for -polarization. The detailed procedure to solve (7) and (8) using the method of 
moments is described in [16]. 

To summarize the proposed method, we first compute the equivalent magnetic current over 
the cavity's aperture by applying the finite element method to (2). Then, we compute the 
surface magnetic field induced on the region of the cavity's aperture by applying the 
method of moments to (7) and (8), or (9) and (10) when there is no coating. Finally, we 
calculate the radiated field using (5) and (6). 


III. Numerical Results 

In this section we present some numerical results to demonstrate the validity and capability 
of the formulation described above. Without loss of generality, the same microstrip patch 
antenna, depicted in Fig. 2, is used in all examples to follow. This antenna consists of a 
rectangular conducting patch fed with a probe and housed in a larger dielectric filled 
rectangular cavity. The first resonance mode at f = 3.3 GHz radiates a field whose £-plane 
is perpendicular to the shorter side of the patch and whose H-plane is perpendicular to the 
longer side of the patch. 


The first set of results is computed for the microstrip patch antenna on a finite circular 
cylinder sitting on a conducting plate, illustrated in Fig. 3. The radiation patterns are 
calculated for two different positions of the antenna whose longer side is parallel to the axis 
of the cylinder. These patterns are shown in Fig. 4 and compared with measured data 
provided by Mission Research Corporation, Dayton, OH. As can be seen, the agreement is 
surprisingly good within the first 30-dB range, demonstrating the validity of the method. 
The disagreement in the deep shadow region is due to the field diffracted by the sides of the 
plate perpendicular to the cylinder and also by the ends of the cylinder, which are ignored 
in our simplified approximate calculation. The corresponding results when the antenna is 
rotated 90 degrees are shown in Fig. 5. Note that in this case, the predicted sidelobes are 
significantly higher that those in Fig. 4 because of stronger surface waves propagating 
around the cylinder. 

The proposed method can, of course, handle multiple patch antennas or arrays, such as the 
one sketched in Fig. 6. The patches and dielectric substrate have the same parameters as the 
one in Fig. 2. The edge-to-edge distance between two adjacent patches is 1.5 cm and the 
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array is centered at <p = 225°. Figure 7 shows its radiation patterns with and without the 
plate. The cavity size for the case of Fig. 7(a) is 12 cm in the 0 -direction and 6 cm in the 
z -direction, and for the case of Fig. 7(b) is 15 cm in the 0-direction and 5 cm in the z- 
direction. The asymmetric pattern in Fig. 7(b) for the case without the plate is due to the 
asymmetric placement of the feed. 

To demonstrate the capability of the method as well as to show the effect of dielectric 
coating, we consider a single patch antenna on a circular cylinder having a diameter of 
30.48 cm. The antenna is the same as depicted in Fig. 2 and the cylinder is coated with a 
layer of material having a thickness of 0.546 cm, relative permittivity e r = 2.2 and relative 
permeability /i r = 1.2 for the lossless case, and e r = 2.2 - 7 1.0 and jU r = 1.2 - y' 0.5 for the 
lossy case. The results are presented in Fig. 8 for the cylinder without coating, with a 
lossless coating, and with a lossy coating. As can be seen, when the longer side of the 
patch is parallel to the axis of the cylinder, it does not excite a surface wave propagating 
around the cylinder. Therefore, the coating does not have a significant effect on the 
radiation pattern because it is not thick enough to support a surface wave. Flowever, when 
the shorter side of the patch is parallel to the axis of the cylinder, it excites a surface wave 
around the cylinder, which is further enhanced by the presence of a lossless coating, 
raising the level of radiation in the shadow region. When the coating is lossy, the surface 
wave is attenuated, resulting in a lower level of radiation in the shadow region. 

To show the effect of a finite ground plane, we consider the microstrip patch antenna of 
Fig. 2 placed at the center of a flat conducting plate whose width is w and thickness is 0.5 
cm. The radiation pattern in the plane perpendicular to the plate is shown in Fig. 9(a) when 
the longer side of the patch is parallel to the axis of the plate and in Fig. 9(b) when the 
shorter side of the patch is parallel to the axis of the plate. The proposed technique provides 
a simple way to study the finite ground plane effect on microstrip antenna radiation patterns 
[17]. 

Finally, we consider a waveguide-fed trihedral, depicted in Fig. 10, to show the 
applicability of the proposed method to any cavity-backed conformal antennas. The 
waveguide is 0.9 inch wide and 0.4 inch thick and is at the center of one of the planes. The 
radiation patterns at f - 9.0 GHz in both the E- and /7-planes are given in Fig. 1 1 along 
with the measured data. Good agreement is observed over a wide range of angles. The 
disagreement in the shadow region and in the region close to 90° in Fig. 1 1 (a) is mainly 
due to the third plane, which is ignored in our simplified approximate calculation. Since the 
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waveguide radiator is not a highly resonant structure, we note that the external structures 
close to the waveguide can have nonnegligible effect on the field distribution at the 
waveguide opening. As a result, the proposed method can be less accurate than in the case 
of highly resonant microstrip patch antennas. 


IV. Conclusion 

In this paper, we described a method to calculate the radiation patterns of cavity-backed 
microstrip patch antennas on a cylindrical body of arbitrary cross section that may be coated 
with dielectrics. In this method, we first employed the finite element method to characterize 
the microstrip patch antennas and then applied the reciprocity theorem, in conjunction with 
a two-dimensional method of moments, to calculate the radiated fields. Several numerical 
examples were presented and compared with measured data to demonstrate the validity of 
the method. Numerical results were also given to show the effects of patch orientation, 
dielectric coating, and finiteness of a ground plane on the radiation patterns. The presented 
method can be used to investigate the effects of external structures (such as the shape and 
coating) on the radiation patterns of microstrip antennas and to provide reference data for 
validating more general methods for calculating the radiation patterns of microstrip antennas 
on general three-dimensional bodies. 
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Figure 1 . Cavity-backed microstrip patch antennas on the surface of a cylindrical 
body with an arbitrary cross section. 
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Figure 2. Configuration of a cavity-backed microstrip patch antenna, (a) Top 
view, (b) Side view. 
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Figure 3. A microstrip patch antenna on a circular cylinder attached to a plate, (a) 
Three-dimensional view, (b) Cross-sectional view. 
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Figure 4. Radiation patterns of a microstrip patch antenna on a circular cylinder 
attached to a plate, with the longer side of the patch parallel to the axis of the 
cylinder, (a) For antenna 1. (b) For antenna 2. 
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(a) 



(b) 


Figure 5. Radiation patterns of a microstrip patch antenna on a circular cylinder 
attached to a plate, with the shorter side of the patch parallel to the axis of the 
cylinder, (a) For antenna 1 . (b) For antenna 2. 
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(a) 



(b) 

Figure 7. Radiation patterns of three microstrip patch antennas depicted in Fig. 6. 
(a) With the longer side of the patches parallel to the axis of the cylinder, (b) With 
the shorter side of the patches parallel to the axis of the cylinder. 
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(a) 



Angle (degrees) 


(b) 


Figure 8. Radiation patterns of a microstrip patch antenna on a dielectric coated 
circular cylinder, (a) With the longer side of the patch parallel to the axis of the 
cylinder, (b) With the shorter side of the patch parallel to the axis of the cylinder. 
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(a) 



Angle (degrees) 


(b) 


Figure 9. Radiation patterns of a microstrip patch antenna on a flat conducting 
plate having width w. (a) With the longer side of the patch parallel to the axis of the 
plate, (b) With the shorter side of the patch parallel to the axis of the plate. 
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Figure 10. Configuration of a waveguide-fed trihedral. The waveguide is 0.9 
inch wide and 0.4 inch thick and located at the center of one of the planes. 
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APPENDIX 2 


Hybrid FEM/SBR Method to Compute the Radiation Pattern from a 
Microstrip Patch Antenna in a Complex Geometry 

A. D. Greenwood, S. S. Ni, J. M. Jin, and S. W. Lee 
Center for Computational Electromagnetics 
Department of Electrical and Computer Engineering 
University of Illinois at Urbana-Champaign 
Urbana,, Illinois 61801-2991 

Key terms: Finite element method, shooting and bouncing ray method, 
microstrip patch antenna, radiation. 

Abstract 

A previously introduced technique to compute the radiation of a microstrip patch 
antenna on a cylindrical body of arbitrary cross-section is extended to include patch 
antennas mounted on complex, three-dimensional bodies. The technique involves using 
the finite element method (FEM) to characterize the patch antenna, and then comput- 
ing the radiation pattern using the reciprocity theorem and the shooting and bouncing 
ray (SBR) method with edge diffraction. The validity of the calculation and the need 
for an extension to three dimensions is shown by comparison with measured data for a 
waveguide-fed trihedral. Other calculations show the effects of the external geometry 
on the radiation pattern of a microstrip patch antenna. 

1 Introduction 

Recently, a technique for calculat ing the radiation pattern of a microstrip patch antenna 
on a cylindrical body was developed [1]. The technique is a hybrid method using 
the finite element method (FEM) and the two-dimensional (2-D) method of moments 
(MoM). The FEM is used to analyze the patch antenna with its aperture closed by a 
perfect conductor, resulting in an equivalent magnetic current density, M. The 2-D 


25 



MoM is then used to find the radiation due to M in the presence of the cylindrical 
body. While this technique is simple and works well, its application is limited to patch 
antennas on 2-D bodies such as cylinders. Although the technique can be extended to 
3-D bodies by employing a 3-D MoM, the computational cost is prohibitive. Thus, an 
alternate method of extension is needed. 

In this paper, a technique is presented which replaces the 2-D MoM with the high 
frequency shooting and bouncing ray (SBR) method combined with edge diffraction. 
A microstrip patch antenna is now considered on a 3-D complex body. As in the 
original technique, the FEM is first used to analyze the antenna with its aperture 
closed by perfect conductor, yielding an equivalent magnetic current density, M, over 
the aperture. The SBR method combined with edge diffraction is then used to find the 
radiation due to M in the presence of the 3-D complex body. The proposed method 
can also be regarded as the extension of the hybrid FEM/SBR scattering method [2] 
to radiation problems. 

This paper first discusses the theoretical development of the method. It then 
presents some numerical results, showing the utility of the method. Finally, it compares 
results generated by the extended technique to measured results for a waveguide-fed 
trihedral, showing the validity of the method and the need for a 3-D extension. 

2 Theoretical Development 

Consider a microstrip patch antenna mounted on a 3-D complex structure. Since a 
microstrip patch antenna is a highly resonant structure, its current distribution and 
input impedance at resonance are determined mainly by its internal structures such 
as the shape of the patch and the thicknesses and the dielectric constants of the sub- 
strate and the superstrate. Provided that the patch antenna is placed on a locally flat 
surface, the shape of the host 3-D body has little effect on its current distribution and 
input impedance. Therefore, the patch antenna can be characterized approximately by 
assuming that it is placed on an infinitely large ground plane and using the FEM as 
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described in [3]. The FEM analysis provides the solution to the equivalent magnetic 
current density, M, over the aperture of the patch antenna. 

Once the equivalent magnetic current density has been calculated, the radiation 
pattern is obtained by computing the field radiated by M in the presence of the host 
3-D body. An accurate method to compute the radiated field is to use the 3-D MoM. 
However, the computational cost is prohibitive, especially for a large 3-D body. A more 
efficient approach is to use a high-frequency approximation such as the SBR method 
[4, 5]. 

There are two approaches to using the SBR method to compute the radiated field. 
One approach is to first compute the radiated field over a small hemisphere covering 
the patch antenna. This field is then converted into many rays, which shoot along 
the radial directions. Each ray is traced as it bounces around the 3-D body, and and 
the bounces are governed by geometrical optics (GO). At the last hit point, or at 
each and every hit point, a physical optics type integration is performed to determine 
the ray contribution to the radiated field. The final result is the summation of the 
contributions from all the rays. This approach, called the forward calculation, has the 
advantage of simultaneously computing the radiated field in all directions. However, 
to obtain accurate results, the field on the hemisphere surrounding the patch antenna 
must be divided into many rays, and in tracing each ray, its divergence factor must be 
calculated and tracked. 

The second approach is to employ the reciprocity theorem. In this approach, called 
the backward calculation, an infinitesimal dipole is placed at the observation point, 
and the SBR method is applied to the dipole radiation. Because the dipole is located 
far from the structure containing the patch antenna, it produces a plane wave incident 
upon that structure. The SBR method computes the fields radiated by the dipole using 
GO and accounting for the phase shifts and the polarization properties of the fields. 
Because the incident field is due to a plane wave, the ray divergence factor is set to 
unity. Denoting the field radiated by the patch antenna as Emspa, the field radiated 
by the dipole as Hd, and the current density on the dipole as J, the fields Emspa and 
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Hd are then related by the reciprocity theorem, 

JJJ v Emspa ' J dV = - jf H d • M dS 

where V denotes the volume of the dipole and 5 denotes the aperture of the patch 
antenna. The integral on the left hand side of the equation evaluates to the component 
of Emspa in the direction of J times the strength of the dipole, II. The quantity 
Hd on the right hand side of the equation is also proportional to 77, which is then 
canceled from the expression. Two orthogonal components of Emspa are found by 
considering both horizontal and vertical polarizations for the incident plane wave in 
the SBR calculations, corresponding to different physical orientations of the dipole 
antenna. 

When using the SBR method to find Hd on the surface of the patch antenna, there 
may be some directions for which the incident wave is blocked. The SBR calculation 
will then give Hd= 0, an erroneous result. To find an accurate value for Hd, edge 
diffraction must be included with the SBR calculation. 

The geometry for the diffraction computation is shown in Figure 1. An arbitrary 
edge is subdivided into short segments, and the diffracted electric field from a short edge 
segment is computed using the “incremental length diffraction coefficients” (ILDC) 
introduced by Mitzner [6, 7]. The diffracted field from the entire edge is then computed 
from the sum of the contributions of the segments, and the diffracted magnetic field 
is found from the electric field using Maxwell’s equations. Given a short edge segment 
of length I , the diffracted field at a distance $ from the segment is given in terms of 
ILDC by [6, 7] 

Ej{P) 1 __i dee d e<t> E' e (Q) 

E$(P ) J s [d# d H \ E'^Q) 



28 



in which d uv is given by 


d — — e 3 ke (' s -' s ')'i 

U UV — A C 

4ir 


(u ■ i) (t ■ v) D s 4- it ■ (t x s)(t x s') ■ v F^j 


sine 


k(s - s') ■ t 


( 2 ) 


where sinc(a;) = sinz/a;, and D s and Dh are diffraction coefficients. The coefficients 
D s and D h. are defined in terms of the Keller-Ufimtsev diffraction coefficients as 


D s = 2 e J>/4 
D h = 2e J>/4 


2nk 

sin /3' 0 sin /?o 

2irk 


sin 3' 0 sin (3q 


D 

D 


dif 


dif 
h * 


(3a) 

(3b) 


Denoting the interior wedge angle as a = (2 - n) it, Df { and Df l{ are given by [8, 9] 


Df = r (F>i - D 2 ) 

2nJ2nks\n P' Q sin j3 o 

(4a) 

e ~ J,r / 4 

Ff = T (D 1 + D 2 ) 

(4b) 

2nJ2itks\\\ /3 f 0 sin (3q 


where 


D\ = cot 
Z ?2 = cot 


TT + (0 - 1 />') 


2n 

\tt + (V> + V’O 


2rc 


F[&Z,< 7 + (i/> - ^')] + cot 
F[/cZ,£r + (^ + V>')] + cot 


7T - (V> - Ip') 


2n 

f 7T - (V> + ^') 


2ft 


F[ArFfir (i/> - V>')] (5a) 

F[fcF^ _ (^ + V’')]. (5b) 


In the above, the function F[-] involves a Fresnel integral and is defined by 


/ oo . , 

e-*' dt, 


(6) 


and the other parameters are given by 


L = s sin j3’ 0 sin /3 0 and g ± (X) - 1 + cos (X - 2nnN ± ) (7) 
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where N ± are the integers which most closely satisfy 


2n7riV + - A" = +7T and 2 nnN - X = -n. (8) 

It should be noted that the expressions above compute only the first order edge 
diffracted field. Thus, the incident field used in the calculations is that resulting from 
direct illumination of the edge. Multiply bounced rays are not included in the edge 
diffraction computation, nor are fields diffracted by one edge and then subsequently 
diffracted by a second edge. In most cases, a first order edge diffraction computation is 
adequate for determining the dominant effects on the radiation pattern of a microstrip 
patch antenna. 

3 Numerical Results 

To illustrate the usefulness of the proposed method, numerical results are generated for 
a patch antenna in two different geometries. The patch antenna under consideration 
for these tests is shown in Figure 2. The antenna resonates at 3.3 GHz, and, according 
to the proposed method, it is characterized by an equivalent magnetic current density 
found using the FEM. 

The first example involves placing the patch antenna on a finite, circular ground 
plane. In Figure 3, the resulting computations of the H-plane and E-plane radiation 
patterns are compared to the patterns generated by the same antenna on an infinite 
ground plane. As can be seen, significant features of the antenna pattern are missed by 
assuming an infinite ground plane to make calculations. Edge diffraction is primarily 
responsible for the differences between the finite ground plane examples and the infinite 
ground plane results. 

For the second example, the patch antenna is placed on a stair-step geometry shown 
in Figures 4a and b. In Figures 4c and d, the computed radiation patterns (H-plane 
and E-plane) are compared to those found by assuming an infinite ground plane to 
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make the calculations. Like the finite, circular ground plane case, significant features 
are missed by assuming an infinite ground plane, and edge diffraction is primarily 
responsible for these features. 

To show the validity of the method, numerical computations are compared to mea- 
sured data for a waveguide-fed trihedral. The geometry for this comparison is shown 
in Figure 5a where the waveguide is 0.9 inch wide by 0.4 inch thick, and the operating 
frequency is 9.0 GHz. Under these conditions, an incident TEio wave propagates in 
the waveguide, and the aperture is characterized by an equivalent magnetic current 
density, similar to the characterization of a patch antenna by an equivalent current 
density. The calculated radiation patterns (H-plane and E-plane) are compared to the 
measured patterns in Figures 5b and c. The calculated results agree with the measured 
results very well. 

Although it is a 3-D problem, this waveguide-fed trihedral configuration was consid- 
ered previously using the 2-D MoM [1]. The MoM result agrees well with the measured 
result for the H-plane pattern when the azimuth angle is between -90° and 40°, and 
it agrees well for the E-plane pattern when the elevation angle is between -40° and 
90°. However, it misses the oscillations in the radiated power outside of these regions. 
The oscillations that are missed in the 2-D MoM computation are correctly predicted 
by the SBR/edge diffraction computation and are caused by edge diffraction in the 
third dimension. Thus, this problem not only validates the method, but underscores 
the need to extend the original 2-D technique to consider 3-D problems. 


4 Conclusion 

The technique for computing the radiation pattern from a microstrip patch antenna 
on a cylindrical body [1] was extended to compute the radiation from a patch antenna 
in a 3-D complex geometry. The method first involved using the FEM to analyze the 
antenna. The result was an equivalent magnetic current density on the aperture of 
the antenna. Then, the SBR method was used to calculate the radiated field, and for 
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this a plane wave was launched toward the antenna by an infinitesimal dipole. The 
SBR method was used to find the resulting magnetic field on the aperture of the patch 
antenna, and the radiated field was calculated using the reciprocity theorem. To insure 
accuracy, the field diffracted by edges was computed using ILDC and added to the SBR 
result. Numerical results showed that the structures surrounding the patch antenna 
have dominant effects on the radiation pattern, and assuming an infinite ground plane 
to make calculations misses these effects. Comparison of numerical results generated 
by the method to measured data showed good agreement. Finally, comparison of nu- 
merical results from the method presented here to those from the original 2-D MoM 
technique showed that pattern features missed by the 2-D technique are correctly pre- 
dicted by this technique. 
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Figure 1: Geometry of edge diffraction. Note that ip' and xp are measured from an 
illuminated face of the wedge, so 0 < xp f < n and 0 < xp < n7r where a = (2 — n)n . 

Figure 2: Microstrip patch antenna considered in computational tests. 

Figure 3: Radiation pattern of a microstrip patch antenna on a finite, circular ground 
plane. 

Figure 4: Radiation pattern of a microstrip patch antenna on a stair-step geometry. 

Figure 5: Comparison of calculated data with measured data for the radiation pattern 
of a waveguide-fed trihedral at 9.0 GHz. 
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SCATTERING FROM A CYLINDRICALLY CONFORMAL 
SLOTTED- WAVEGUIDE ARRAY ANTENNA 

Guo-Xin Fan and Jian-Ming Jin 
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ABSTRACT 

A numerical method is developed to investigate electromagnetic scattering 
by a cylindrically conformal waveguide-fed slot array. The problem is first for- 
mulated in terms of integral equations using the equivalence principle. The 
integral equations are then solved using the method of moments (MoM) in con- 
junction with global sinusoidal basis functions and Galerkin’s testing procedure. 
The MoM solution requires the evaluation of the generalized admittance ma- 
trices involving various dyadic Green's functions. The slow convergence of the 
series associated with the summation of waveguide modes is accelerated using 
the Kummer transformation and the slow convergence of the series associated 
with the summation of the exterior modes is avoided by using the asymptotic 
solutions with proper treatment for singular integrals. The evaluation of the 
excitation vector and scattered field is also accelerated using Watson's transfor- 
mation and asymptotic solutions. Numerical results are presented to illustrate 
the scattering characteristics of the cylindrically conformal waveguide-fed slot 
arrays, such as the effects of curvature, slot thickness, and waveguide termina- 
tion on the radar cross section of the arrays. 
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I. INTRODUCTION 


Slotted-waveguide array antennas are widely used on modern airborne radars 
because of their large power handling capability, high efficiency, light weight, 
compact structure, and low side lobes. However, being highly efficient radiating 
structures, by reciprocity theorem these antennas are also efficient scatterers, 
which contribute significantly to the overall radar cross section (RCS) of the host 
vehicle. Usual RCS reduction techniques cannot be applied to these antennas 
without degrading their performance. Bandpass radomes, made of frequency 
selective surfaces (FSS), can reflect nearly all the incident energy outside the 
working frequency band and allow the incident energy within the working fre- 
quency band pass through to reach the antennas. However, it is at the working 
frequency band that the slot arrays have a significantly higher RCS. Therefore, 
it is necessary to predict the RCS of the slot arrays for applications such as 
target identification, electromagnetic compatibility, and stealth technology. 

Because of a large number of slots and the mutual coupling between the 
slots through the waveguides and exterior space, a full-wave analysis of slotted 
waveguide arrays is very difficult. In the past, most work was focused on the 
radiation analysis of a single slot [1]- [ 1 0] , one-dimensional arrays and small two- 
dimensional arrays [11]-[17], all in a planar surface. In particular, Stevenson 
[1] developed what is now considered the classical theory for a single slot and 
Oliner [2] presented a variational solution of the problem. Khac and Carson [3] 
employed the method of moments (MoM) to seek a numerical solution to the 
slot field using pulse basis functions and point match technique. The efficiency 
of the MoM solution was improved by Lyon and Sangster [5] and Stern and 
Elliott [6] by using global and piecewise sinusoidal basis functions and Galerkin’s 
technique. The MoM was also applied to tilted slots [9], dielectric-covered slots 
[8], [10], edge slots [18], and slots in a sectoral waveguide [19]. Recently, Fan [20] 
analyzed cylindrically conformal slotted-waveguide array antennas with a curved 
waveguide as the feeding guide and sectoral guides as the radiating guides. 
Whereas the literature for the radiation analysis of the slot arrays is abundant, 
little work was done on the scattering from the planar slot antennas [21]-[23], let 
alone the curved ones. Josefsson [21] analyzed scattering by a single slot in an 
infinitely long waveguide and Chen and Jin [23] developed the MoM and finite 
element method (FEM) solutions for the slots in the waveguides terminated 
with arbitrary loads. No work was found on the analysis of the scattering 
characteristics of slotted waveguide arrays on a curved surface. In this paper, 
we present such an analysis. 
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The key problem in the analysis of slotted waveguide arrays is to solve for 
the slot aperture fields. Among various numerical methods, the MoM is most 
efficient and accurate for this purpose because only a small number of basis 
functions are needed to represent the aperture fields. However, the major diffi- 
culties in the MoM are (i) its formulation requires the dyadic Green’s functions 
for the waveguide, the exterior space, and the slot if it has a finite thickness, 
and (ii) the evaluation of its matrix involves highly singular integrals and very 
slowly converging series. In this paper, we address all these problems. To be 
more specific, we first formulate the integral equations for the problem and 
apply the MoM with global basis functions and Galerkin’s testing procedure. 
We then describe in detail the evaluation of the MoM matrix involving various 
dyadic Green’s functions. The slow convergence of the series associated with 
the summation of waveguide modes is accelerated using the Kummer transfor- 
mation and the slow convergence of the series associated with the summation 
of the exterior modes is avoided by using the asymptotic solutions developed by 
Boersma and Lee [24] and Bird [25]. The evaluation of the excitation vector and 
scattered field is also accelerated using Watson's transformation and asymptotic 
solutions. Finally, we present some numerical results to illustrate the scattering 
characteristics of the cylindrical lv conformal slotted waveguide arrays, such as 
the effects of curvature, slot thickness, waveguide termination, and frequency 
on the RCS of the arrays. 


H. THEORY 

In this section, we describe in detail the formulation of the problem and its 
solution by the MoM. Particular attention is given to the computation of the 
elements of the generalized admittance matrices. 

A. Integral Equations and MoM Solution 

Consider the problem of electromagnetic wave scattering by a waveguide-fed 
slot array on a cylindrical surface whose radius is p 3 , as depicted in Fig. 1. The 
cross section of each waveguide is an annular sector with inner radius p \ , outer 
radius p 2 , and subtended angle A <j>, and the thickness of the curved walls of the 
waveguides is t = p 3 — p 2 - Each waveguide is terminated at z = z\ and zn with 
an arbitrary load. All radiating slots are longitudinal slots cut in the outer wall 
of the waveguides, having the same width 2w and different length and offset 
with respect to the center-line of the waveguides. 

For the ith slot, its inner and outer apertures, S( and Sf , divide the space 
into three regions: the waveguide region (region a), the region outside the cylin- 
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der (region 6), and the slot region (region c), as illustrated in Fig. 2. In ac- 
cordance with the equivalence principle, the fields in the three regions can be 
decoupled by covering the apertures S( and Sf with a perfectly conducting 
surface and introducing equivalent magnetic currents above and below the per- 
fectly conducting surface. Denoting the equivalent magnetic current below S( 
as Mf and that below Sf as Mf , because of the continuity of the tangential 
electric fields across the apertures, the equivalent magnetic current above S{ is 
— Mf and that above Sf is — Mf . Therefore, the field in region a is due to Mf, 
the field in region b is due to —Mf and the incident field, and the field in region 
c is due to — M f and M f . 

By enforcing the continuity of the tangential magnetic fields across Sj and 
Sf , we obtain the integral equations satisfied by the equivalent magnetic cur- 
rents, Mf and Mf , as 


^H?,(Mj) + H e Ti (Mf)-m,(Mf) = 0 

j 

res' 

(1) 


r € Sf 

(2) 


j 


where the subscript r denotes the tangential component, the summation in (1) 
is carried out for all inner apertures in the same waveguide, and the summation 
in (2) is carried out for all outer apertures. The magnetic field is related to the 
surface magnetic current by 

H°(r) = JJ 5“(r,r') • M(r')dS' (3) 

5 " 

where a == a for region a. b for region 6, and c for region c, respectively, and 
G’ a (r,r') denotes the magnetic-source magnetic-field dyadic Green's function in 
the corresponding region. Finally, H prt is the field due to the incident field in 
the presence of the conducting cylinder without slots. 

To seek a numerical solution of (1) and (2), we first expand each equivalent 
magnetic current using the global sinusoidal basis functions 


M j = zlkfj = :£ V ti siaa vt 

r 6 S' 

(4) 

9 = 1 






M f = £Mf = i£v* sin 

resf 

(5) 

q=l 




where a q j = q^/Lj with Lj being the length of the j th slot, and choose the 
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testing or weighting functions as 

Wp f = i Wp,- = £ sin ctpiZ res' (6) 

Wp,- = = isinop,^ r £ S? . (7) 

Applying Galerkin's procedure, the integral equations can be converted into the 
matrix equation given by 

■ [y«(o) + y;j(c,s!-s<)} [-y#(c,s/ ; s/)] 1 r [vr] j non 

. [-Y#(C'Sf ;5/)j [y«(6) + Yji(c,S?;Sf)}\ \ [l# J \ [I pi ] J 

( 8 ) 

where [V^ (a)] and [Y^ (6)] are the generalized admittance matrix for regions 
a and 6, respectively. |y p ^ ^c, S- is the generalized admittance matrix for 

region c, in which S{ is the aperture on which the field point is located and 
Sj is the aperture on which the source point is located. Finally, [I p i] is the 
generalized excitation vector. 

We note that for slots of zero thickness, the edge condition is often employed 
to represent the transverse distribution of the equivalent magnetic current [6], 
[10], [14], which makes the evaluation of the matrix elements difficult. For slots 
of finite thickness as is the case for practical applications, the singularity is much 
weaker than that for the zero thickness case. As a result, the equivalent magnetic 
current can be expanded in terms of the global sinusoidal basis functions with a 
uniform transverse distribution. Such an expansion has been successfully used 
in the past for the analysis of rectangular waveguide slot radiation problems [5], 
[8], [9], 

In the remainder of this section, we describe the evaluation of the elements 
of each generalized admittance matrix and excitation vector, necessary for a 
numerical solution of (8). 


B. Generalized Admittance Matrix Elements for Region a 

The admittance matrix elements due to the internal coupling in the sectoral 
waveguide (region a) are given by 


’”< o)= IJJI G - 


(r, r') sin a P i£ sin a q j£ dS dS ' 


where G?. is the if-component of the magnetic-source magnetic-field dyadic 
Green s function for the sectoral waveguide, which is derived in Appendix A 
and given by 


G a zz = jun Y B(m..n) C^< 

m = Q 
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Also, in the above, u = mn/Ao, CA = cos i/y>, C ' v > = cos i/y?' , /i m „ = ,/A.' 2 — A,- nn , 
Amn are the roots of B[,(X mn po ) = 0, A m „ = 1 - R< ln R> ln £- jh "' n{ ‘ 3 ~ z '\ and 
R< n and R> n are the reflection coefficients for mode (?n, ?/) at c = zi and 
c = ^ 2 , respectively. For shorted terminals R mn — — 1, and for matched termi- 
nals R mn = 0. 

The integrals in (9) can be evaluated analytically through the transformation 
of variables, resulting in 

fL, ru rLj rll 

/ / / 6 f ?,(r, r') sin cr r ,£ sin a qj ^ clj/d^'drjd^ 

3 J — w J 0 J — w 


Y£(a) = 


CO CO 


= -jui€ Y2 H fi (m, n)I'J(m) [/! J 0 (m. n) + 7‘ 2 m (»n,n)j (14) 

m = 0 n = 0 
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m ± 0 (16) 


and rf,- denotes the offset of the ith slot from the center of the waveguide. Also, 
/i’o(m,n) = C(m.n) [l + (-i)W] [l - ( -i)P e -J*-^.] + a~ P7 

®pi ^mn 

(17) 

/* q ( TH , n ) = C7 ( 771 , T3 ) hm r.[(s°‘ — ) — { -L. — i-j )/2] 
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with 


_\ _ mn "Pf MQ\ 

C(m,n)_ 2 _ J, 2 a 2 _ 1.2 

and coi and ro j denote the center of the ith and jth slots, respectively. Finally, 
/^ m (m,n) can be written in three parts as 

n ) + Ci( m - n ) (20) 


( m , n ) = — C(m, n)/J^ n e _jAm "(^' 0,+io ^ (£.+£j )/ 2 2 -d 
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While the summation in (14) can be evaluated without any difficulty when 
i jL j, its evaluation when i = j involves a slowly converging series given by 

CO 

S(m) = B(m } n)/lo — Si(m) + 5o(m) (25) 

n = 0 

where 

CO 
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n = 0 
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Slim) = , _ 7 . (27) 

n = 0 C V 

To show this, let us examine the asymptotic behavior of the series when n — > oo. 
Using the asymptotic expression of the Bessel functions for large arguments, we 
obtain the approximate solution to the eigenvalue equation B t l/ {X m np 2 ) — 0 as 

A mn ~ — — (28) 

P 2 - P\ 
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Clearly, the asymptotic form of the terms in S\(m) is proportional to n -3 and 
that in 52 (m) is proportional to n” 2 . To accelerate the convergence of 52(m), 
we apply the Kummer transform [26] to (27), resulting in 


S 2 (m) = S p<l (k 2 - al,)LA £ 

lm=o L ~r 


B(m,n) f p 2 — p\ V B(rn. cx 

~ h m n V n ) n2 + x 
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where 
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( 31 ) 


(32) 


The first summation in (31) converges rapidly and the second summation in 
(31) can be written in a closed form as 


00 1 
«- + x “ 
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+ 7T COth TTx] 
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In passing, we note that when pi and pn oc, the expressions for sectoral 
waveguides reduce to those for the corresponding rectangular waveguides with 
the following substitutions 


P2_ _ 

a 


(34) 

V 

777 7T 


x n 

/ m7T \ - 

/ 71 7T \ 2 

(35) 

Knn ~ 

(-) 

+ (t) 

B(m,n) ~ 

2 

1 

(36) 

1 + <^0m 

k 2 ab 


where a = p\ A<$ or p 2 ^<^ and b = p 2 — p\ are the dimension of the rectangular 
waveguide. 


C. Generalized Admittance Matrix Elements for Region b 

The admittance matrix elements due to the external coupling in region b are 
given by 

Vp q { b)= fj JJ G b ..(r,r')sma F i£ sin a<jj£' dS dS 1 (37) 

s' s' 

• J 

where G b zz is the ii-component of the magnetic-source magnetic-field dyadic 
Green's function for the conducting cylinder. The rigorous expressions of the 
dyadic Green’s function involve infinite series and infinite integrations of Han- 
kers functions [27], [28]. Although these expressions can be converted into the 
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forms suitable to numerical calculation [28], it is still very time consuming to 
compute the matrix elements, especially when the cylinder is large. Therefore, 
it is necessary to seek the fast converging asymptotic solutions. Five different 
approximate asymptotic solutions have been developed in the past and these 
are denoted as PINY [29], UI [30], OSU [31], B-L [24], and TSB [25] solutions, 
and each has its accuracy and range of validity. A comparative study [28] shows 
that the B-L and TSB solutions offer the best overall accuracy, and these two 
solutions are complementary to each other. More specifically, the B-L solution 
is more accurate when the field point is near the source or in the paraxial region 
while the TSB solution is better when the field point is far from the source and 
off axis. 

The B-L solution [24] is given by 


6 ';. ~ Go(0 jw(4) [cos 2 0 + <7 ( 1 — (j) (2 — 3 cos 2 <?)] 

'll 2 


<iQv'(0) 


11 /II 2 187 

- cos- 0 + <1 ( — + - tan- 0 - — cos’ 6 


(38) 


and the TSB solution [25] is given by 

G'j. ~ Go(t) |f(j) [cos 2 6 + q(l - v)(2 - 3 cos 2 0)] - <} v(d) (^sin 2 0 - 


, ,11 17 . , 

+t ' lW f 60 “ 36 Sin " 


+ v 2 (/?) ( sin 2 6 + 


40 ) + 5 0 


( 0 ) 


(39) 


where t = sf(z - ;')*' - (<p — Tp ' )-' , 0 = arctan[(r - r')/(y? - v?')], 9 = j/kt, 
0 = (k cos 4 $/‘2R~) l f 3 t , and v, V\ , v 2 , and Vq 1 * are the surface Fock functions, 
all defined in [24], [25], and Appendix C. 

When i = j, G b zz has a t~ 3 singularity as t — > 0, and a regularization is 
needed to evaluate the integrals in (37). To do this, we first rewrite (38) as [25], 
[30] 

G b z = G° zz + G z: (40) 


where 


G° zs 

G zz 


Go[t ) |cos 2 § + q (1 - q) (2 - 3 cos 2 0 )j 
G 0 (0| [v(4)-l] [cos 2 0 + g ( 1 — g) (2 — 3 cos 2 0)] 
-q0v' (/?) 


11 ,= . /II 2 o- 187 , 

-cos-fl + 7 (- + -tan-tf-— cos -0 


]}• 


(41) 


(42) 
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It can be shown that G?. is the ff-component of the dyadic Green's function 
for an infinite conducting ground plane, given by 






e -jk t 


27 rt 


(43) 


whereas G zz can be considered as the perturbation due to the curvature of the 
cylinder, with a singularity of t ~ 3y/2 as t — > 0. As a result, the self-admittance 
matrix element can be written as 


v;*(6) = Y™(b) + y"(4). 


(44) 


Using integration of part and transformation of variables, V p ° 9 “(6) can be written 
as [8] 

l2 i . / i \p + <? fL* flw e -jkt 

= l <«> 

where t = y/u- + v 2 and 

( T-3. [ Q >”^' 2 ” q 9<) sin °V' U “ - a pi) sin Q r< u ] 


F{u,v)= { 


P 5* 9 

(2u- - 1 >)[(I» - «)(i' 2 - o*,)cosQp,u+ ^-(P 2 +a 2 ,)sina r ,u 

p = q. 


(46) 

The integral in (45) becomes regular when evaluated in the polar coordinates. 

To evaluate the perturbation term in (44), we first let rj = Rp, if = Rp * , 
where /? = p 3 , and write it as 

Yp' q {b) = f f f f G lz (Z- ? ,i]~ rf) sina pi £ sin a q i? dr) 1 d£ dr] d£. 

J 0 J — \jj J 0 J —\u 

(47) 

Introducing the transformations u = £ — t j = tj — 7?', td = £ 4- £' — and 
t/ = rj 4- r/, and observing that G* 2 (u,r) is an even function of u and i>, we 
obtain 

?”(6)= [l + (-l) p+? ] £' r( 2 w-v)G„(u,v)F(v)dvdu (48) 


where 


F(u) = 


q - — a * 

p» <f * 


(- 1 ) p+9 a p, sin a q i u - a,, sin a pi u j 


(Li — u) COS QpiU + ~ sin Q f ,;U 


P 5* 9 
P = <7- 


(49) 
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With another transformation u — tcosO^ and v — £sin0, (48) becomes 


r 1 ( ft* o rLyf cos 6 r2w / sin 6 \ 

,»=[> + (-!)'«]{/ I + JJ o } 

■(2w - / sin 0) G 2Z (t,6) F(t } 9) t d$ dt 


where #0 = tan - 1 (2w/L t ) and G* ; (U#)has the same form as (42). The integrals 
in (50) become regular by letting t = r 2 , which can be evaluated by expanding 
the integrand in terms of the power of r and using Gaussian quadrature. The 
Fock function can evaluated using its small-argument expansions. 

When i ^ j, the expression for the admittance matrix elements can be 
written as 


/*L, r to fL } r W 

YpqW = / / / / G tl (: - z',R(p- >p')) sin a p it sin a qi £' dt)' d£' dr] d£ 

J 0 J — to J Q J — to 


using the coordinate transformation 


= (C-C) + (ro i -Zo J ) + (I,-L J )/2 

R(ip — = fa — rf) + R(tp 0i — ip Q j) 

where ^ 0 « and y> 0j are the coordinates of the center of the ith and jth slot, 
respectively. 

In practical slot arrays, the length difference between the slots is small. In 
this case, we can use the average slot length of the two slots to approximate 
their length, that, is, 

Li zz Lj zz (Li 4 - Lj)/2 = L. (53) 

Using the transformations and u' = £ + - L, and the midpoint 

integration, we obtain 



f [G zz {u 
Jo 


) + (-l) p+ *G j; (-u)] F(u)du 


(54) 


where F(u ) is given by (49) by replacing L, with L, and G :z (u) is given by (38) 
and (39) with the transformations : — = u + (r 0l — r 0 j) and R(ip - p') = 

R{<Poi ~ Foj)- 

In passing, we note that in this work, the TSB asymptotic solution is used 
for two distant slots whose centers have the same 2 -coordinate, whereas the 
B-L asymptotic solution is employed for all other cases. We also note that the 
asymptotic solution of the surface fields on a cylinder can be interpreted as the 
creeping waves propagating from the source to the field point along the surface 
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of cylinder. In general, there are at least two creeping waves. However, for an 
open cylindrical surface, only the shortest one has to be considered. 


D. Generalized Admittance Matrix Elements for Region c 

It can be shown (Appendix B) that the magnetic field in the cavity (region c 
with a sectoral cross-section), formed by the slot with its two apertures covered 
by a perfectly conducting surface, due to the surface magnetic current z6(p — 
p') sin a P ,'£ is given by 


0 - 

H z {p,p') = - — R(p,p') sinap.S 

juip 


where 0- = |i- 2 — o ", | and 


R(P, P ') = 


C B(p,,p) B{p 3 ,p') 
CB(p 3 ,p) B(p 3 ,p’) 


P< p' 
P> P 1 


with 


(55) 

(56) 

(57) 

(58) 

(59) 

0P2 ’ ^ ~ dp3 (60) 

Using (55), we obtain the expressions for the generalized admittance matrix 
elements for region c as 


C = (l/3)[B(p,,p , )B'(p 3 ,p')- B(p 3 ,p')B'(p2,p')Y 
B(p\p) = 


J' o (3p')Y o (0p) - Y£{0p')M0p) k > o r , 
l' o (0p' ) K o (0p) - A'o (Bp 1 ) Io(3p) k< api . 


It is easv to see that 


B'{pi,p 3 ) = B'(p 3 ,p 3 ) = 0 
B{p2, Pi) = -J— , B(p 3 , p 3 ) = 




r<? ' 


pq ' 





B'(p 3l pi) u 1 

LiW0 

-1 

JWH 

3 Pi B'(p2,p 3 ) 

Liwd 

1 


dp 3 B'(p 3 ,p2) 

L t wd 

-B(p 2 ,p 3 ) 


dij dpq 


dij dpq 


(61) 

(62) 

(63) 

(64) 


J jojp B’(p2,p 3 ) 

When p 2 and p 3 — ► oo, the expressions above reduce to those for a rectangular 
cavitv as 


Y& (c,S!-,Sj) = Y*i{c, Sf ; Sf ) = ^cot 3t S tj S Fq 
Y‘ J q (c, 5/ ; 5/ ) = Y;j (c, Sf ; Sj ) = — ^4- S ti 6 


j sin Qt 




(65) 

( 66 ) 
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for k > a pi, and 


Ypq (CtS, 1 ;Sj) = Ypq(c, Sf ; Sf ) 
Ypq ( c i Sf \Sj) = Yp* {c,S? ; Sj ) 

for k < a pt -. 


—r coth 0t SijSpq 

JUJ^I 

Liiv(3 -1 £ £ 

jw/i sinh 0t ij pq 


(67) 

(138) 


E. Excitation Vector 

Consider an incident plane wave of arbitrary polarization and incident angle, 
whose electric field is given by 


E inc = (6 inc cos 0 + p inc sin v)e~> k r (69) 

where v is the polarization angle, (# inc ,^ mc ) denote the incident angles, and 
(r mc , 0 inc , p inc ) are the unit vectors. In the computation of the excitation 
vector, the primary field, H pri in (2), is the sum of the incident field and the 
field scattered by the cylinder in the absence of all the slots. This field is given 
by 


H prt = —Yo sin v sin 0” lc rJ kzcose ' nc M (kpz sin 0 inc , p — <p inc ) 

P-P 3 


where 


M(x,p) = -f-Y — 

V YJ j TTX t-il+di 


2 j n COS Tip 


Jnx ^ 1+don H p w 

From (70) we obtain the elements of the excitation vector as 


(70) 


(71) 


2 P * 


= // 


W r ' sin apii dS 

P=fi 3 


= sin ti> sin 6 inc e J * Zo ' cos s '" e f p (6 ,nc ) 


■M(kp 3 sin 6 ine ,ip 0 i — <p inc ,&<f) 


[72) 


where 


fp(0) = 


~(jb cos6 ) 


i -i 


— cos ( 


kL, 


'>0) 


Ap) = — Y 

it rr ^ ' 


f cos np 


j„ x ^ ol+6on 


Sam(n Ay?) 


p = odd 


j sin ( * cos 0) p = even 


[73) 


(74) 
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in which = u'/pa, and Sam(x) is a sampling function defined as 


Sam(x) = ( . 

f sir 

For narrow slots, we can use the midpoint integration to find 


x = 0 
sin x/x i /0. 


10 


Ipi = 4Yo—^—— sin 4' sinO inc ei kz °- cos ' ,nc f p (6 inc )M(kp 3 sm6 inc ,s oi - <p inc ). 
pn 

(76) 

When x becomes large, Al{x,<p) for the lit region can be evaluated more effi- 
ciently using the expression [27] 


A/ (x, ip) = 2e^ xcostp 


(77) 


which can be considered as the geometrical optics approximation. Note that 
(77) has the same accuracy as Gorianiov’ more complicated expression [32]. For 
the neighborhood of the shadow boundary and the shadow region, M(x, <p) can 
be evaluated using the expression obtained using the Watson transformation 
and Fock theory [27], 




o-j*( *-i) 


(>-§)(?) 


1 / 3 ' 



(78) 


where p(x) is the far-zone Fock function, whose definition is given in [27], [33], 
[34] and Appendix C. 


F. Scattered Field and RCS 

Once the equivalent magnetic currents on the slot apertures are found, we 
can compute the scattered field, which can be considered as the radiated field 
due to the equivalent magnetic currents. Using the Fourier transform and far- 
field approximation [35], it can be shown that the field radiated by Mf is 


E e = 0 (79) 

_ kr Nf /S 

E v = sin e M (kp 3 sin 6, <poi - SP, As?) T (9) (80) 

2n ‘ r p?i P 

where 7‘, = r — zo,cos0. For a narrow slot on a large cylinder, we can again 
use the approximation described in preceding subsection to simplify the far-field 
calculation. 
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The RCS of the slot array is defined as 


2 


<7 = lim 4 t rr 2 

r — too 


E 5 CO (r) 

E ,nc (r) 


where E 3ca is the total scattered field from the slot array. 


(81) 


III. NUMERICAL RESULTS 

The formulation described above has been implemented in a computer pro- 
gram and a variety of numerical results have been obtained to study the effects 
of various factors on RCS. Because of limited space, we present only some typical 
results. 

For all the examples considered, the cross-section of sectoral waveguides is 
given by \{p\ + P2)A<^> = 21 mm and P 2 — pi — 10.16 mm. The distance 
between the center-lines of the adjacent waveguides is 22 mm, measured on the 
surface of the cylinder. The adjacent slots on the same waveguide are spaced 
A^o/2 apart in the r direction, where X g o is the waveguide wavelength at the 
working frequency. The shorted planes are placed A ff o/4 away from the first 
and last slots, respectively. In order to avoid the numerical overflow in the 
computation of the generalized admittance elements for the shorted waveguides 
at the working frequency, the position of the shorted planes is displaced by Ar, 
which can be considered as an acceptable manufacturing tolerance. All other 
parameters, unless otherwise specified, are given by: 2 w = 1.6 mm, Li = 1.6 
mm, Si = 1.5 mm, A t = 0.8 mm, Ar = 0.01 mm, ip = 90°, and the working 
frequency /o = 9.1 GHz. The normal at the center of the array is 0 = 90° and 
0 = 0° in the spherical coordinate system. 

The computer program was first validated by comparing its solution to those 
obtained by MoM and FEM [23] for a single slot on a planar surface. Good 
agreement was observed between the three solutions. 

To demonstrate the effect of the curvature of the host cylinder on the RCS 
of a conformal slot array, Figs. 3 to 6 show the RCS for pi = oo, pi = 500 cm, 
pi = 200 cm, pi = 100 cm, and pi = 50 cm in both E- and H-planes. As can be 
seen, the cylinder’s curvature has a significant effect on the patterns and values 
of the RCS in the E-plane. When p x becomes large, the RCS approaches that 
for a planar array. When pi decreases, the grating lobes begin to disappear 
and the maximum RCS decreases. In contrast, the pattern of the RCS in the 
H-plane remains the same and only the value of the RCS is decreased as pi 
decreases. 
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To further demonstrate the effect of the cylinder’s curvature on the RCS, 
Fig. 7 presents the space distribution of the RCS for p\ = oo, p\ = ‘200 cm, 
pi — 100 cm, and p\ — 30 cm. As can be seen clearly, as p\ decreases the 
grating lobes in the E-plane is dilated where the pattern in the H-plane remains 
the same. 

To show the effect of the waveguide terminals on the RCS, Fig. 8 gives the 
RCS of a 8 x 8 conformal slot array (p\ = 100 cm) when the waveguide terminals 
are matched and shorted. It can be seen that when the plane wave is incident 
normally on the cylinder ( 0 = 90°), the RCS is very small for the array with 
shorted terminals. This is because at the working frequency, the wave entered 
the waveguide is reflected by the shorted terminals and the reflected wave cancels 
the incident wave at the slot. Thus, the total electric field or the total equivalent 
magnetic current is zero at the aperture of the slot. An alternative explanation 
is that with the shorted terminals, the equivalent impedance looking at the 
aperture of the slots is zero. Therefore, when the shorted planes are displaced, 
the equivalent impedance will be changed and so is the RCS. This is clearly 
demonstrated in Fig. 9. 

Figure 10 shows the effect of the slot thickness on the RCS. As can be seen, 
at the working frequency, the slot thickness has no effect on the RCS when 
the waveguides are matched and when the waveguides are shorted, it has a 
significant effect. This effect can also be explained using the concept of the 
equivalent impedance. 

Finally, Fig. 11 shows the frequency responses of a 8 x 8 array (pi = 100 cm) 
with matched and shorted terminals, respectively. It is very interesting to note 
that the RCS of the slot array is almost identical no matter the waveguides are 
matched or shorted when the frequency is not close to the working frequency. 
At the working frequency, the RCS is substantially different. This implies that 
when the frequency of the incident wave is not close to the working frequency of 
the slot array, the energy entering the waveguides is very trivial and the main 
contribution to the RCS is the wave scattered by the slots directly. Therefore, 
in this case the dominant factor is the geometry of the slots, instead of the 
structures behind the slots. Figure 12 shows the space distribution of the RCS 
of a 16 x 16 array (pi = 100 cm) at three different frequencies. 

IV. CONCLUSION 

In this paper, a method of moments (MoM) solution was developed for elec- 
tromagnetic scattering by a cylindrically conformal waveguide-fed slot array. 
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The equivalence principle was first employed to formulate the problem in terms 
of integral equations. Various dyadic Green’s functions required in the formu- 
lation were derived and the integral equations were solved using the MoM in 
conjunction with global sinusoidal basis functions and Galerkin’s testing proce- 
dure. The MoM solution requires the evaluation of the generalized admittance 
matrices involving some slowly-converging series. The slow convergence of the 
series associated with the summation of waveguide modes was accelerated using 
the Kummer transformation and the slow convergence of the series associated 
with the summation of the exterior modes was avoided by using the asymp- 
totic solutions with proper treatment for singular integrals. The evaluation of 
the excitation vector and scattered field was also accelerated using Watson's 
transformation and asymptotic solutions. Numerical results were presented to 
illustrate the scattering characteristics of the cylindrically conformal waveguide- 
fed slot arrays, such as the effects of curvature, slot thickness, and waveguide 
termination on the radar cross section of the arrays. It was observed that (i) 
as the curvature of the host cylinder increases, the grating lobes in the E-plane 
dilate and eventually disappear and the maximum RCS is reduced; (ii) the inter- 
nal structure of the waveguides such as the waveguide terminals has negligible 
effect on the RCS when the frequency of the incident wave is not close to the 
working frequency of the slot array; (iii) at the working frequency the RCS of 
the slot array with a shorted terminals at normal incidence changes significantly 
with the position of the shorted terminals; and (iv) at the working frequency 
the slot thickness has no effect on the RCS of the array with matched terminals 
whereas for the array with shorted terminals its effect is noticeable. The numer- 
ical method developed are useful in applications such as target identification, 
electromagnetic compatibility, EMP penetration, and stealth technology. 



56 



and boundary conditions 


hg m = 0, hxVxg n 


= 0 . 


(84) 


In this work, we require only the rr-component, g m:z , which satisfies the 
scalar wave equation 


1 d f d \ Id 2 8 2 a 

1 r, 1 * Too 4 “ a a + 


L pdp\dpj p 2 d<p 2 , d: 2 

and boundary conditions 

dg mz 


g m!! = - A{p - (85) 


dp 


0, 


dg„ 


d<p 


y3 = 0, 


= 0. 


Using the orthogonality relation [20] 

rP2 

/ pB^iX mn P) B v (A mn •P)dp = ^nn' N \ mn , 

J P i 


( 86 ) 


(87) 


we can show that 

8(p-p' 


6(<p-<p') = y y - V- —4 — 
in, in 1 + ^Orn A<?A Am „ 

( 88 ) 


m = U n =0 


Assuming that 


9m 


» = EE 


1 


' “ 1 -f iom A6AA mri 

m=0n=l) 

and substituting into (85), we obtain 

d 2 


^(A m nP)^(A nm ^)C , v ^VZ 7nn (r 1 r / ) (89) 




(90) 


Enforcing the boundary conditions at z = z\ and ro, which are specified in 
terms of the mode reflection coefficients, R^ n and i?> n , respectively, we obtain 


Zmn(^z) = 


1 




2jh mn A m n 

__ p> j/i r nn(: + ^') i d< d> (i 3 - Cl ) 

- rL mn c c ' iL mn iL mn c c 

=(™J 


91) 


The required rr-component G H can then be obtained from 

(92) 
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with the aid of the differentiation formula for generalized functions 


d 7 

d. 


-jh mn \ 


: “ 2 '' = !*-*'! =2jh mn S(z - z') + 

ozoz 


APPENDIX B 

FIELD IN THE SLOT-CAVITY 

The magnetic field in the slot-cavity, produced by the surface magnetic cur- 
rent zS(p — p')sina p i£, is related to the magnetic vector potential by 


jj _ k 2 / 1 d 2 \ , 

H z — . (id" ) Amz 

jufi\ k-d£ 2 J 

where A mz satisfies the wave equation 

(V 2 + k-)Amz = -6{p - p') sin a pi ^ } r £ V' c 

and boundary conditions 


dA m z 


dp 


= 0, 


dAr , 


<9<p 


= 0, 


dm: 


ly3 = ± A;? 


£=o,z,, 


= 0. 


Assuming A m2 (p,<p,£) = i?(p, p') sin a P ,C we obtain 




df? 

dp 


= 0 . 


(94) 


(95) 


(96) 


(97) 


(98) 


lp = P2 ,P3 

Equation (97) is recognized as the inhomogeneous Bessel's equation of order 
zero, whose solution is given by (55)-(58). 


APPENDIX C 
FOCK FUNCTIONS 

In this appendix, we give the definition and related formulas of the Fock 
functions used in this paper. 

C.l The Surface Fock Functions 
i) Definition : For a complex z and a real x, 

e j(n + l/2)jr/2 ;r n + l/2 


Vn(z) = 


2r(n+ 1/2) 


[ t n ^fie-^dt. n = 0 , 1,2 

Jr i “ 4(0 


(99) 
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e J71 


t s l U x \ - c x / w w> 
° ( ' ~ 20F y Pl [u/ 3 (0 

v(*) = vo(*)i 

v'(x) = [i/(x) — t>j (a:)] /2x, 


e~ Jxt dt, 


where 


w 2 (t) = -7= / exp (*r - z 3 /3) dr. 

>A dr 2 


The integration contour Ti (or To) goes from 00 to 0 along the line specified by 
ang t = — 27t/ 3 (or angr = 27r/3) and from 0 to 00 along the real axis. 

ii) Small- Argument Asymptotic Expansion : For x < 1, 

10 

»n(*) = 1 + (-1)" £ (104) 

m = 1 

10 

^(x) = v / 7re _J ' r/4 x 1/ ' 2 1 + ^ c m (-jx) 3m/2 , (105) 

rn= 1 

where the coefficients 6 rn (0), 6 m ( 1 ) , 6 m (2), and c m are given [25]. 

iii) Residue Series Asymptotic Expansion: For x > 1, 

j(n-l/2)r/2 x n + l/2 

= arm 4 1/2) (,06) 


ti^fx) = \/ffe J ' 3ff/4 x 1/2 ^(1 +Ji-C)^7-3- 


where = |o^ n |exp(-j7r/3), and — are the zeros of derivative of Airy’s 
function, Ai'(x). 


C.2 The Ear-Zone Fock Functions 
i) Definition: 


s(*) = 4 =/ ~n7\ c ~ jxt dt " 

V* Jr l 


ii) Formulas: 


g(x) = 2e"-' rJ/3 (x < -3) 




-3 < x < —0.5) 


= 1.39937+ E %(Kx) m (—0.5 < x < 0.5) 

/ — ' 77) I 
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(0.5 < x < 4.0) 


20 


= V 

a mAl( — °m) 

_ j g 35e -(0 8823-jO 5094 ) J :-jxV3 (l > 4.0) 

where k — exp(— jbw/6), d m , A m and s m are given in Table I. 


( 109 ) 


Table 1 


in 

dm 

^1 m 

Sm 

\ 

0.7473831 

—0.346754038 ID — 13 -j0.2020905526D-12 

-0.8253328597D+01 

jO. 12 193 14385 D+02 

o 

-0.6862081 

— 0.1 469465991 D— 08 -j0.2834631597D-08 

-0.6157249604D+01 

j0.9659289583D+01 

3 

-2.9495325 

-0.3540987938D+00 j0.1274055399D+00 

— 0.4741573876D+00 

j0.2711S64904D+0l 

4 

-3.4827075 

0.2544 176553D— 05 -j0.1407419787D-05 

— 0. 447412691 7D +01 

j0.7632294399D+01 

5 

8.9378967 

— 0.23 123408 15D— 03 J0.4465889350D-03 

-0.3001 183618D+01 

j0.5858670248D+01 

6 

56.1946214 

0.1332955205D— 01 -j0.2046671398D-01 

— 0.1673939722D+01 

j0.4236757852D+01 

7 


0. 1 1 8094480 1D+01 j0.1164437839D+01 

0.572814962 1D+00 

jO. 1 234 148084 D +01 

8 


0.7689249066D+00 -j0.1285238403D+01 

0.1413116779D+01 

— j0.3362910S94D+00 

9 


— 0.2105540865D+00 j0.1312640573D-01 

0.2125953705D+01 

— j0.2263304268D+01 

LO 


0.1622S92013D-02 -j0.1090163259D-02 

0.311763791 1D+01 

— j0.5238642524D+01 
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FIGURE CAPTIONS 


Figure 1. Conformal slotted- waveguide array. 

Figure 2. Division of three regions and surface magnetic currents. 

Figure 3. The effect of the cylinder’s curvature on the RCS in the E-plane for 

a 16 x 16 slot array with matched terminals at / = 9.1 GHz. (a) planar, 

Pi - 500 cm, (b) p x = 200 cm, (c) p x = 100 cm, (d) pi = 50 cm. 

Figure 4. The effect of the cylinder’s curvature on the RCS in the H-plane for 

a 16 x 16 slot array with matched terminals at / = 9.1 GHz. (a) planar, 

p l = 500 cm, (b) p\ — 200 cm, (c) p\ — 100 cm, (d) p\ — 50 cm. 

Figure 5. The effect of the cylinder’s curvature on the RCS in the E-plane for a 

16 x 16 slot array with shorted terminals at / = 9.1 GHz. (a) planar, 

pi = 500 cm, (b) pi = 200 cm, (c) pi = 100 cm, (d) pi = 50 cm. 

Figure 6. The effect of the cylinder’s curvature on the RCS in the H-plane for a 

16 x 16 slot array with shorted terminals at / = 9.1 GHz. (a) planar, 

pi = 500 cm, (b) pi = 200 cm, (c) pi = 100 cm, (d) pi = 50 cm. 

Figure 7. RCS (dBsw) of a 16 x 16 slot array with shorted terminals at / = 9.1 
GHz. (a) planar; (b) pi = 200 cm; (c) pi = 100 cm; (d) pi = 30 cm. 

Figure 8. The effect of the terminals on the RCS of a 8 x 8 conformal slot array 

(pi = 100 cm). • • ■ matched terminals, shorted terminals, (a) E-plane 

pattern; (b) H-plane pattern. 

Figure 9. The effect of the position of the shorted plane on the RCS of a 8 x 8 

conformal slot array (pi = 100 cm). A: — 0.001cm, * - • A z ~ 0.1cm, 

A z = 1cm. (a) E-plane pattern; (b) H-plane pattern. 

Figure 10. The effect of the slot thickness on the RCS of a 8 x 8 conformal slot 

array (pi = 100 cm) in the E-plane. At = 1 mm, • ■ At = 0.4 mm, 

At = 0.1 mm. (a) matched terminals; (b) shorted terminals. 

Figure 11. RCS of a 8 x 8 conformal slot array (pi = 50 cm) as a function of 
frequency, (a) 0 inc = 90°, ^ nc = 0°; (b) 9 inc = 90°, <p ,nc = 45°. 

Figure 12. RCS (dBsw) of a 16 x 16 conformal slot array (pi = 100 cm) at three 
frequencies, (a) / = 8 GHz; (b) / = 9.1 GHz; (c) / = 12 GHz. 
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Figure 1 : Conformal slotted-waveguide array. 

Figure 2: Division of three regions and surface magnetic currents. 
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Figure 3: The effect of the cylinder’s curvature on the RCS in the E-plane 

for a 16 x 16 slot array with matched terminals at / = 9.1 GHz. (a) 

planar, Pi = 500 cm, (b) pi - 200 cm, (c) p x = 100 cm, (d) p x = 50 

cm. 
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(a) 



( b ) 




Figure 4: The effect of the cylinder's curvature on the RCS in the H-plane 

for a 16 x 16 slot array with matched terminals at / = 9.1 GHz. (a) 

planar, pi = 500 cm, (b) p x = 200 cm, (c) p x = 100 cm, (d) p x = 50 

cm. 
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(c) 



(d) 

The effect of the cylinder's curvature on the RCS in the E-plane 

x 16 slot array with shorted terminals at / = 9.1 GHz. (a) 

pi = 500 cm, (b) pi ~ 200 cm, (c) p l = 100 cm. (d) p Y = 50 
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(b) 



(0 



(d) 

Figure 6: The effect of the cylinder's curvature on the RCS in the H-plane 

for a 16 x 16 slot array with shorted terminals at f — 9.1 GHz. (a) 

planar, p x = 500 cm, (b) p x = 200 cm, (c) p x = 100 cm, (d) p x = 50 

cm. 
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(a) 



(b) 


Figure 8: The effect of the terminals on the RCS of a 8 x 8 conformal slot 

array (pi = 100 cm). • • • matched terminals, shorted terminals, (a) 

E-plane pattern; (b) H-plane pattern. 


71 




60 



(a) 



(b) 

Figure 9: The effect of the position of the shorted plane on the RCS of a 8 x 8 

conformal slot array (pi = 100 cm). A z = 0.001cm, • * * Az = 0.1cm, 

Az = 1cm. (a) E-plane pattern; (b) H-plane pattern. 
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(b) 


Figure 10: The effect of the slot thickness on the RCS of a 8 x 8 conformal 

slot array (p\ — 100 cm) in the E-plane. At = 1 mm, • * • At = 0.4 mm, 

At = 0.1 mm. (a) matched terminals; (b) shorted terminals. 
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Figure 11: RCS 
of frequency, (a 




of a 8 X 8 conformal slot array (p i = 50 cm) as a function 
6 inc = 90°, 9 mc = 0°; (b) 6 inc = 90° t = 45°. 
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SCATTERING FROM A LARGE PLANAR SLOTTED 
WAVEGUIDE ARRAY ANTENNA 


Guo-Xin Fan and Jian-Ming Jin 
Electromagnetics Laboratory 
Department of Electrical and Computer Engineering 
University of Illinois at Urbana-Champaign 
Urbana, Illinois 61801-2991 


ABSTRACT 

The problem of scattering from a large planar slotted waveguide array 
antenna is investigated. The structural scattering and slot scattering are 
solved independently through the application of the equivalence principle. 
The former is calculated by a high-frequency method, such as the physical 
optics (PO) and the shooting and bouncing ray (SBR) method combined 
with edge diffraction. The latter is solved using the method of moments 
(MoM) in conjunction with global sinusoidal basis functions and Galerkin’s 
testing technique. Special attention is paid to the evaluation of the ad- 
mittance matrix elements, which involve the internal and external mutual 
coupling. Numerical results are presented to demonstrate the effects of a 
variety of factors on the radar cross section (RCS) of a slot array, and these 
include the number of basis functions in the MoM solution, the external mu- 
tual coupling between the slots, the waveguide terminations, and the host 
object. Also given are the RCS spare distribution and the time responses of 
slot arrays to illustrate their unique scattering characteristics. The RCS of 
a complicated real slot array is given finally to demonstrate the capability 
of the numerical method as well as the distinct features of its RCS pattern. 

I. INTRODUCTION 

Large planar slotted waveguide array antennas are the most common an- 
tennas used on modern airborne radars. Because of their unique structures, 
their radar signature has a number of very unique characteristics. The study 
of these characteristics is important for applications such as target identi- 
fication, electromagnetic compatibility, and stealth technology. Recently, 


76 



we presented a method of moments (MoM) solution to scattering from a 
cylindrically conformal slotted waveguide array antenna [1]. The formula- 
tion is complicated because of the complex nature of the problem. Although 
the method can also be applied to planar slotted waveguide array antennas 
by treating these antennas as the limiting case of a cylindrically conformal 
array antenna with its radius approaching infinity, the solution is unneces- 
sarily complicated and time consuming. Because of the practical importance 
of planar slot arrays, in this paper we present a specific MoM solution to 
this problem based on much simplified formulas, with an emphasis on the 
analysis of large arrays. 

The scattering from a slotted waveguide array antenna consists of two 
parts: the slot scattering and the structural scattering. The former is con- 
tributed by the slots, and the latter is contributed by the surface of the 
antenna and the host object. In practical applications, slot arrays are usu- 
ally placed in the region where only the direct incident field reaches. Be- 
cause of this, the slot scattered field and the structural scattered field can 
be solved separately through the application of the equivalence principle. In 
this paper, we employ a high-frequency method, such as the physical optics 
(PO) and the shooting and bouncing ray (SBR) method combined with edge 
diffraction, to take into account the structural scattering and the MoM to 
characterize the slot scattering. 

The paper is organized as follows. Section II presents the formulation for 
the analysis of planar slot arrays, and the emphasis is given to the calculation 
of all kinds of the generalized admittance matrix elements. Section III gives 
several typical numerical examples to illustrate the scattering properties of 
planar slot arrays, such as the effect of the aperture modes and external 
mutual coupling, the RCS space and probability distribution, the effect of the 
host object or platform, and the time response of a slot array. A complicated 
real slot array is considered finally to demonstrate the capability of the 
numerical method as well as the distinct features of its RCS pattern. 

II. THEORY 

Consider the planar slotted waveguide array antenna depicted in Fig. 1. 
The waveguides have a width a and height 6, and their upper walls have a 
thickness t . Each waveguide may have a different length and may be divided 
into several sections. Each section is terminated with an arbitrary load. All 
the radiating slots are longitudinal and have the same width 2w; however, 
they may have a different length and a different offset with respect to the 
center-line of the waveguides. 

In accordance with the equivalence principle, the field outside the an- 
tenna can be decoupled from those inside the antenna by covering the slots 
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Figure 1: A planar slotted waveguide array antenna. 


with a perfect conductor and introducing an equivalent magnetic current 
over the apertures of the slots. Therefore, the total scattered field can be 
decomposed into two parts. One is the field scattered by the surface of the 
antenna with all the slots covered by a perfect conductor and this is of- 
ten referred to the structural scattered field (or structural mode), denoted 
as E 5 *™. The other part is the field radiated by the equivalent magnetic 
currents and this is often referred to as the slot scattered field (or antenna 
mode), denoted as E s/ot . The total scattered field can then be written as 

E 5ca (r) = E stru (r) + E slot {r). (1) 

When the surface of the antenna is illuminated only by the incident 
wave and when there is no multiple reflection between the antenna and its 
host object, the structural scattered field and the slot scattered field can be 
solved for separately. Since there are many efficient high-frequency methods 
available to calculate the structural scattered field, we focus on the solution 
of the slot scattered field. 

A. Integral Equations and MoM Solution 

Consider the ith slot. Its lower and upper apertures, denoted as S{ and 
Sf respectively, divide the space into three regions: the waveguide region 
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(denoted as region a), the external region (denoted as region 6), and the 
slot region (denoted as region c) since the slot has a finite thickness. By 
imposing the continuity of the tangential magnetic fields on S( and Sf, we 
obtain the integral equations satisfied by the equivalent magnetic currents, 
and M f , as 

£ H“t(M|) + H^-(Mf) - H‘,(Mf) = 0 res! (2) 

3 

£ H 6 ri (M/) + W ri (Mf) - = H?7 r € Sf (3) 

3 


where Mf denotes the equivalent magnetic current on S- and M f denotes 
the equivalent magnetic current on Sf , H a denotes the magnetic field in 
region a, H pri denotes the primary field which is the sum of the incident 
field and the field scattered by the antenna without the slots, and finally, the 
subscript r denotes the tangential component of the associated field. Note 
that the summation in (2) is carried out for all inner apertures in the same 
waveguide, and the summation in (3) is carried out for all outer apertures. 

For narrow slots, we can choose the global sinusoidal basis functions to 
expand each equivalent magnetic current 


TV 


ij 


M}' 1 = zMj'* = i 


Y V 1 : 

g3 

<?=1 


IJ 


sin aqj€ 


re Sj 


IJ 


(4) 


I jr 

where a qj = qn /Lj with Lj being the length of the j th slot, Nj denotes 
the number of basis functions, and £ varies from 0 at the one end of the 
slot to Lj at the other end. By employing Galerkin’s technique, the integral 
equations can be converted into the matrix equation given by 


•[y«(«) + y«(c,s/;Sj)] 



( K\ } 



Y^(b) + Y^(c,Sf;Sf)}_ 

\ [V3] J 

1 [/*■] I 


(5) 


In the equation above, the generalized admittance matrix elements and ex- 
citation vector elements have the same physical meaning as those of the 
conformal slot array [1]. 

We note that for slots of zero thickness, the edge condition is often em- 
ployed to represent the transverse variation of the equivalent magnetic cur- 
rent [2]-[4], which makes the evaluation of the matrix elements difficult. For 
slots of finite thickness as is the case for practical applications, the singularity 
is much weaker than for the zero thickness case. As a result, the equivalent 
magnetic current can be expanded in terms of the global sinusoidal basis 
functions with a uniform transverse distribution. Such an expansion has 
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been successfully used in the past for the analysis of rectangular waveguide 
slot radiation problems [5]- [8]. 


B. Generalized Admittance Matrix Elements for Region a 

The admittance matrix elements due to the internal coupling in the rect- 
angular waveguide (region a) are given by 


w = //// r') sin a p i£ sin dS dS f 


(6) 


si s] 


where G a zz is the ii-component of the magnetic-source magnetic-field dyadic 
Green’s function for the rectangular waveguide, which is given by 


OO OO s 

G a zz = jut n ) C x C x 'C y C y ,U{z - z’) 

m — Dn — fl ' 

1-2— j g — jhmn + ) 


m=0n=0 

Ai 


h-mn A 

-/?> n e- jUm " z V' lmn < 2+2 '> + (1 - A mn )e j/l ”*"l 2 - 2 'l] (7) 

In this expression, 

2 2 1 


J B(ra, n) = 


(8) 


1 *T 1 T ^On 

and C x — cosvx, C x > = cosz/x', C y = cos^t/, C y / = cos/xt/, = m7r/a, /x — 
W&. A^ n = ^ 2 + /z 2 , /i mn = v/^ - A^ n , A mn = (**“*•), 

and f?< n and Rf nn are the reflection coefficients for mode (m, n) at the two 
ends of the waveguide, located at z = z\ and z = 22? respectively. 

The integrals in (6) can be evaluated analytically through the change of 
variables, resulting in 

rLi rw rLj rw 

Ypq (a) = / / / / G a zz (r, r') sin a p ,£sin otqj^'drj'd^'drjd^ 

J 0 J—wJO J —w 

OO OO 

= Z Z n ) p x ( m ) [n j o ( m - «) + ( m > n )] ( 9 ) 


m=0 n=0 


where 
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( 10 ) 


(k 2 — a 2 )L 

Izo( m , n) = C(m, n) [l + (-l) p+? ] 1 - (-1 * 

L J CC. — 


h 2 

pt ,v mn 
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( 11 ) 

Zoi <; ZQj (12) 
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with 


C(m, n) = 




a 


pi 




qj 




pi 


h 2 a 2 — h 2 

n mn ** qj n ™ 


(13) 


in which 8 X denotes the offset of the ith slot from the center of the waveguide, 
and zoi and zoj denote the center of the ith and jth slots, respectively. 
Finally, I l z \ m (m, n) are the same as those in [1]. 

For the mutual coupling between different slots (that is, when i ^ j), 
the summation of the double series in (9) can be evaluated directly without 
any difficulty due to the fast evanescence of higher-order waveguide modes. 
For the self coupling element (that is, when i = j), its evaluation involves a 
slowly converging series, given by 


~ (&-<&) Li k 2 - - 2 

5Z fi ( TO - n ) 5 P1 ~~2 _ h 2 = 6 P1 L > 

n=0 a pi n 'mn 


k 2 ab 


* W ^ 


1 


+ fion n 2 + X 2 


(14) 

where x 2 = ( b/n) 2 {a 2 i - k 2 + v 2 ). However, this term can be evaluated 
analytically in a closed form since 


V' 2 1 * , 

?J+<o„»= + ^ = ^ COth ' 1 ' 


(15) 


n=0 


C. Generalized Admittance Matrix Elements for Region b 

The admittance matrix elements due to the external coupling in region 
b are given by 

Y pq ( b ) = fj JJ G h zz (r, r') sin a pi £ sin a qj €' dS dS' (16) 

Sf sj 

where G zz is the zz-component of the magnetic-source magnetic-field dyadic 
Green’s function for the outer space. As usual, when the surface of the 
antenna is large enough and the interference from the surrounding objects 
can be ignored, G zz can be approximated by the corresponding Green’s 
function G zz for an infinite conducting ground plane, 


g: 


Gt = 


k 2 


1 + 


i d 2 


p -jkt 


k 2 dz 2 I 2nt 


(17) 


where t — yj(x — x 1 )^ + (z — z*)*. 

When the source point and the field point are in the same aperture (that 
is, when i = j), G zz has a £“ 3 singularity as t — >• 0, and a regularization is 
needed to evaluate the integrals. However, this singularity can be reduced 
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to the order of t 1 using integration by parts. As a result, Y^(b) can be 
written as 


F(u,v) 


■■ k 2 l + (-l) p + 9 f f $0 r L */ cos 0 m r/2r2w/sm6^ 

pg jio\i 27 rk 2 1 Jo Jo + ] 9o Jo J 

xF(£sin 6,t cos 6) e~^ ki dtdO (18) 

where Oo — tan -1 (2w/L{) and 

f 2(2 w — u) r 2 9 \ ' n2 2 \ • 1 

2 ~ 2 [ a pi\ k ~ a qi) Sin a qi U “ <*qi{ k ~ <*pi ) Sin a pi U \ 

Q pi ® qi 

p^q 

(2 10 - w) (If - u)(k 2 - Op,) cosQpfM + -^-{k 2 + a pt ) sin a p ;u] 

p = q. 

(19) 

In arriving at (18), we have applied a series of variable transformations [7], 
When i / j, the expression for the mutual admittance matrix elements 
can be written as 

rL x rw rLj rw 

Y pq ( b )= / / / G„{x-x',z- z') sin a p ,£ sin a gi ^ dr{ d £' dr) d £ 

JO J—wJ 0 J —u) 

(20) 

with 

2 - 2 ' = (£-£') + ( 2 o. - 2 0 j) + (I,- - Ij)/2 r9 . 

*-*' = (l?-*?0 + (Xoi-*Oi) 1 ’ 

where xo t and xoj are the coordinates of the center of the ith and jth slots, 
respectively. For large practical slot arrays, the length difference between the 
slots is small. Using the average slot length of the two slots to approximate 
their lengths [7], [8], that is, 


Li % Lj % (Li + Lj)/2 = L, 


( 22 ) 


and introducing the transformations u = £ — and v! = £-f — L ) with the 

aid of the midpoint integration, we have 

yjg (b) = 2 w 2 J q L [Go(tt) + (-l) p+9 G 0 (-u)] F(u) du (23) 
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D. Generalized Admittance Matrix Elements for Region c 

The procedure similar to that in [1] can be employed to obtain the ex- 
pressions for the generalized admittance matrix elements for region c. The 
magnetic field in the cavity formed by the slot with its two apertures covered 
by a perfectly conducting sheet, produced by the surface magnetic current 
z8(y — y f ) sin a p if, is related to the magnetic vector potential by 


„ k 2 / 10 2 \ 

— I 1 T i 0 Qr9 1 ^ mz 
k 2 o£ 2 J 

where A mz satisfies the wave equation 

(V 2 + k 2 )Amz = -S{y - y') sin a pi 
and the boundary conditions 


re v c 


dA r 


dx 


= 0, 


dA v 


= 0, 


y=0,f 


Assuming A mz (x, y,£) = Y (y, y') sin a pi f , we have 

d 2 Y 


€=0 ,L, 


= 0. 


dy 2 


+ I3 2 Y = -6(y-y r ), 


dY_ 

dy 


= 0. 


y=o,t 


(26) 

(27) 

(28) 

(29) 

(30) 


where (3 = yjk 2 - a 2 pi for k > a pt -, and (3 = ~ k 2 for k < a pi . The 

solution of (29) and (30) is 


Y{y,y ') = - cos (3y < cos (3(t — y>) //?sin f3t (31) 


where t/< = min(y,y'), t/> = max(y,y'). Evaluating the inner product be- 
tween the aperture magnetic fields and the testing functions, we finally ob- 
tain 

y; 3 q (c, sU sj ) = y ;> ( C , s, f ; s/) = - —cot pt s tJ s pq (32) 

JOJ/J, 

Y£(c , 5/;5/) = y«(c, 5/; 5/) = csc/3t (33) 


E. Excitation Vector 

For an incident plane wave of arbitrary polarization and incident angle, 
whose electric field is given by 

E ,nc = ( 0 mc cos $ + 4> inc sin ^)e- jk '" c - r (34) 
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where 0 is the polarization angle and (0 ,nc , <fi inc ) denote the incident angles, 
we obtain the elements of the excitation vector as 

hi = 4 J Y Tk 2 [(- 1 ) p e Ail - l] smiPsm6 inc smc{k x w) 

^ J/0 {Z0i ~ Li /2)] (3^) 

where y 0 denotes the y coordinate of the upper apertures of the slots. Note 
that for scattering calculation, the spherical coordinate system is defined 
such that the normal at the center of the array is 0 = 90° and <f> = 0°, 
which is obviously different from the rectangular coordinates used for the 
waveguides. 


F. Scattered Field and RCS 

Once the equivalent magnetic currents on the slot apertures are found, we 
can compute the scattered field by slots, which is the radiated field due to the 
equivalent magnetic currents. Meanwhile, for a large array, we can employ 
a high-frequency method, such as PO or SBR, to calculate the structural 
scattered field. The total scattered field is the sum of these two fields. 

The RCS of the slot array antenna is defined as 


a = lim 47rr 2 

r — foo 


E sca (r) 

2 

— 1 1 tti Arrr^ 

E s, °*(r) + E* <ru (r) 

E ,nc (r) 

— 11 111 ‘i/i / 

r— »oo 

E mc (r) 


(36) 


Another parameter of practical importance is the RCS probability of a slot 
array antenna over a space solid angle £7 0 , which can be defined as 


P{cr > a 0 ) = 


jj- JJ U[a(0,(f)) — gq] sin 6 d6 d<f) 

Ho 


(37) 


where <r 0 is a given value of RCS, and U(x) is the unit step function. 


III. NUMERICAL RESULTS 

For all the examples considered except for those mentioned otherwise, 
the cross-section of waveguides is given by a = 21 mm and b = 5.0 mm. 
The distance between the center-lines of the adjacent waveguides is 22 mm. 
The adjacent slots on the same waveguide are spaced X g o/2 apart in the z 
direction, where A p0 is the waveguide wavelength at the working frequency. 
The shorted planes are placed A fl0 /4 away from the first and last slots, re- 
spectively. All other parameters, unless otherwise specified, are given by: 
2 w = 1.6 mm, L t = 16 mm, <5, = 1.5 mm, t = 0.8 mm, ip = 90°, and the 
working frequency /o = 9.1 GHz. 
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To verify the theory above and the computer program, comparisons are 
made between its solution and those obtained by MoM and FEM [9] for a 
single slot and four slots on a single waveguide. Good agreement is obtained 
between the three solutions. In addition, the validation is also confirmed 
by comparisons between the results of a planar slot array and those of a 
conformal slot array on a very large cylinder as shown in [1]. 

We first study the effect of the number of basis functions or aperture 
modes in (4) on the calculation of RCS. The choice of this number depends 
on several factors, such as the slot offset, the height of waveguide, and fre- 
quency. In large practical slot arrays, the slot offset is usually small. Numer- 
ical results show that within a finite bandwidth near the working frequency, 
only three modes are needed for a good accuracy. In many cases, especially 
at a frequency near or below the working frequency, a single-mode approxi- 
mation can yield an acceptable accuracy. Figures 2 and 3 show the effect of 
the number of modes on the RCS of a slot array (excluding the structural 
scattered field) at three frequencies. The use of a minimum number of modes 
makes it possible to analyze very large slot arrays. 

Next, we study the effect of the mutual coupling through the exterior 
space (region b). Unlike the internal mutual coupling that only exists be- 
tween the slots in the same waveguide, the external mutual coupling exists 
between all the slots. Although the fourfold integral in (16) has been re- 
duced to a single integral in (23) for the computation of the external mutual 
coupling, it is still very time consuming to calculate all these matrix ele- 
ments when the slot array is very large. One possible solution is to ignore 
the mutual coupling when the two slots are far apart. Figure 4 shows the 
effect of the external mutual coupling on the RCS patterns (excluding the 
structural scattered field). Three curves are plotted in the figure: one with 
all the mutual coupling, another with no mutual coupling at all, and the 
third one with the coupling with the surrounding 9x5 slots (9 waveguides 
each containing 5 slots). It is apparent that by including the mutual coupling 
between the adjacent slots, one can obtain a sufficient accuracy. 

The effects of the structural and slot scattered fields are illustrated in 
Fig. 5, which shows the RCS pattern of a circular slot array consisting of 
32 waveguides with 32 slots in the central waveguide of the array (hereafter 
referred to as 32 x 32 circular slot array). The radius of the antennas is 
16 in. The pattern of the antenna without the slots (that is, the structural 
scattered field) is also given for comparison. It can be seen that the structural 
scattering forms a strong peak in the normal direction of the plate. The 
major effect of the slot scattered field is the grating lobes. 

Figure 6 gives the space distribution of the RCS of the 32 X 32 circular 
slot array, including both the structural and slot scattered fields. The main 
RCS peaks away from the normal direction are due to the slot scattering. 
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(a) 



(b) 



(c) 

Figure 2: The effect of the number of aperture modes on the RCS in the 

E-plane for an 8 x 8 slot array with shorted terminals. 7 modes; 

5 modes; 3 modes; ••• 1 mode, (a) / = 8 GHz; (b) / = 9.1 GHz; (c) 

/ = 12 GHz. 
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(a) 



(b) 



(c) 

Figure 3: The effect of the number of aperture modes on the RCS in the 

H-plane for an 8 x 8 slot array with shorted terminals. 7 modes; 

5 modes; 3 modes; ••• 1 mode, (a) / = 8 GHz; (b) / = 9.1 GHz; (c) 

/ = 12 GHz. 
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(a) 



(b) 

Figure 4: The effect of the external coupling on the RCS for a 32 x 32 

rectangular slot array with shorted terminals at / = 9.1 GHz. all slots; 

9x5 slots; no external coupling, (a) E-plane; (b) H-plane. 


From the data in Fig. 6, we can calculate the RCS probability using (37) and 
the result is shown in Fig. 7. The space solid angle flo in (37) is specified 
by a cone half-angle, denoted as 9 0 , where 0 o = 0 in the direction normal to 
the surface of the antenna. 

To demonstrate the effect of a slot array antenna on the RCS of the host 
object, we place the 32 x 32 circular slot array antenna in the nose of an 
airplane. The radome covering the antenna is ignored because it is largely 
transparent at the working frequency. The airplane is in the normally flying 
position (horizontal) and the RCS of the entire object is given in Fig. 8. 
The structural scattered field is calculated using the high-frequency code 
XPATCH [10], which is based on the SBR method, and for the calculation, 
the airplane is assumed to be a perfect conductor. We note that although 
the size of the slot array is small compared with that of the airplane, its 
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(a) 



(b) 

Figure 5: RCS of a 32 x 32 circular slot array with matched terminals at 
/ = 9.1 GHz. (a) E-plane; (b) H-plane. 


scattering dominates in some angles corresponding to those of the grating 
lobes. 

The time response of a slot array is studied using a slot array of 8 identical 
waveguides with 2 slots on each waveguide as depicted in Fig. 9. The time 
response, shown in Figs. 10 and 11, is obtained by using the inverse fast 
Fourier transform (FFT) of the frequency response over a bandwidth from 
8 to 12 GHz. In the example shown in Fig. 10, the waveguide terminals 
are matched and the two different intervals between the slots are used, d = 
100 and d = 200 cm, respectively. The first two peaks correspond to the 
scattering directly from the two rows of the slots, respectively, and their 
intervals in time are 2dcos0 mc jc — 3.3 and 6.6 ns, as expected. The third 
peak comes from the first interaction between the slots: the wave enters the 
waveguide from the first (and second) slot and comes out from the second 
(and first) slot. In the example of Fig. 11, the one end of the waveguides 
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Figure 6: The RCS (dBsm) space distribution of a 32 x 32 circular slot array 
with shorted terminals at / = 9.1 GHz. 



Figure 7: The RCS probability distribution of a 32 X 32 circular slot array 
with shorted terminals at / = 9.1 GHz. 








T 









d 



Figure 9: An 8 x 2 slot array. 



Figure 10: The time response of a 8 X 2 slot array with matched terminals 
exited by an incident wave with 6 inc = 60°, (f> tnc = 0°. (a) d = 100 cm; (b) 
d = 200 cm. 

is matched and the other end is shorted at three different positions. The 
second response is due to the field entered the waveguide through the slots, 
reflected by the shorted terminal, and re-radiated from the slots. It is clearly 
seen that the time delay of the second response is proportional to the length 
of the waveguides. 

In the last example, we consider a real slot array antenna constructed in 
our laboratory. The configuration of one quarter of this antenna is shown 
in Fig. 12. Each waveguide is divided into several short sections, each con- 
taining only four slots except for those next to the edge of the antenna. The 
entire antenna consists of 40 waveguides with 42 slots in the central waveg- 
uide of the array (hereafter referred to as 40 x 42 circular slot array). Other 
parameters are: a — 21.3 mm, 6 = 5 mm, 2 w = 3.3 mm, Li = 16.2 mm, 
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Figure 11: The time response of a 8 X 2 slot array with d = \ 3 q/2 and a 
shorted terminal at z x and a matched terminal at z i exited by an incident 
wave with 9 inc = 90°, 4> inc = 0°. (a) z x = 50 cm; (b) z x = 100 cm; (c) 
Z\ = 150 cm. 

Si — 2 mm, t = 0.5 mm, and the working frequency / 0 = 9.9 GHz. The 
calculated RCS patterns are given in Fig. 13, which include the structural 
scattered field. As can be seen, the RCS pattern in the H-plane is substan- 
tially different from that when the waveguides are not divided into several 
sections [see Fig. 5(b)], whereas the pattern in the E-plane remains similar. 
There are several additional grating lobes in the H-plane, which are caused 
by an array whose elements are the short sections. Finally, Fig. 14 shows 
the RCS when this antenna is placed in the nose of the airplane depicted in 
Fig. 8(a). 

IV. CONCLUSION 

This paper investigates the problem of scattering from a large planar 
slotted waveguide array antenna. The total scattering consists of two parts, 
namely the structural scattering and slot scattering, and these two parts are 
solved independently through the application of the equivalence principle. 
The former is calculated by a high-frequency method, such as the PO and 
SBR method. The latter is solved using the MoM in conjunction with global 
sinusoidal basis functions and Galerkin’s testing technique. Special attention 
is paid to the evaluation of the admittance matrix elements, which involve 
the internal and external mutual coupling. Numerical results are presented 
to demonstrate the effects of a variety of factors on the RCS of a slot array, 
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Figure 12: An 40 x 42 circular slotted-waveguide array. 


and these factors include the number of basis functions in the MoM solution, 
the external mutual coupling between the slots, the waveguide terminations, 
and the host object. Also given are the RCS space distribution and proba- 
bility curves and the time responses of slot arrays to illustrate their unique 
scattering characteristics. The RCS of a complicated real slot array is given 
finally to demonstrate the capability of the numerical method as well as the 
distinct features of its RCS pattern. The numerical method developed are 
useful in applications such as target identification, electromagnetic compat- 
ibility, EMP penetration, and stealth technology. 
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(a) 



(*>) 

Figure 13: RCS of a 40 X 42 circular slot array with shorted terminals at 

/ = 9.9 GHz. (a) E-plane; (b) H-plane. 
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